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ABSTRACT 

Three computer simulations of command to line-of-sight 
missile guidance systems were developed within an electronic 
warfare environment to test the feasibility of improved 
guidance techniques while maintaining simplicity for generic 
application among the various missiles of this general 
guidance classification. The first simulation modeled the 
basic guidance scheme; the second introduced a "lead-angle" 
concept; and the third simulation combined these techniques 
for use depending upon whether the target took evasive 
action or employed electronic countermeasures . It was found 
that consideration should be given for use of a "lead-angle" 
variant in conjunction with the basic guidance technique to 
enhance the effective engagement envelope of these missile 
systems against relatively slower-maneuvering targets employing 
only low-duty-cycle denial jamming. 
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I. 



INTRODUCTION 



Command to line-of-sight guidance for several types of 
surface-to-air missiles was developed through computer 
simulation. Consideration was given to electronic counter- 
measures against the guidance methods of this type. 

The simulations are of sufficient complexity to demonstrate 
feasibility while maintaining simplicity for generic application 
among the various missiles of this general guidance classi- 
fication. 
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II. OVER-VIEW OF GUIDANCE TYPES 



Guidance systems can be considered to fall, generally, into 
one of three basic types: line of sight, pursuit, or propor- 

tional navigation. 

The line-of-sight system maintains the missile on the vector 
from the ground-based point of missile control (the tracker) 
to the target aircraft. Variations to this system can provide 
for "lead-angle" intercept through missile command computer 
generation of a synthetic line of sight upon which the missile 
is commanded to fly for target intercept. Also, beam-riding 
missiles fit this category. However, in the beam-riding type of 
line-of-sight guidance the missile avionics are more complicated 
and expensive because the missile normally senses its position 
relative to the beam center and generates its own guidance 
commands internally. This extends the missile's conceivable 
range but usually necessitates some form of terminal homing. 
Basic command to line-of-sight systems can utilize relatively 
cheap missiles by performing guidance computations in the 
tracker unit and not in the missile itself. This limits the 
missile to short-range scenarios and requires that the missile 
possess a greater speed advantage over the target. 

Pursuit guidance systems maintain the missile's velocity 
vector toward the target. Generally, sensors in the missile 
can be assumed to perform this function; therefore, the missile 
can be forgotten by the launcher after the missile has locked 
onto the target. This type of system is capable of medium 
range work against a slowly-maneuvering target. It has 
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relatively simple processing avionics. 

Proportional navigation guidance systems hold a constant 
angle be twe en the mis sile-to- target line of sight and the 
missile axis, thereby, generating a "constant-bearing-decreasing- 
range", collision situation. In other words, the change in 
missile heading is adjusted proportionally to the rate of 
change in t he m issile-to-target line of sight. This type of 
guidance is good for long ranges (and short) against highly- 
maneuvering targets. However, the missile requires either 
seeker and angle-rate sensing, or a command guidance signal 
based upon a collision course. 



III. TYPICAL ENGAGEMENT SEQUENCE 



In order to provide a "vehicle" through which to better 
understand the basic aspects of command to line-of-sight 
guidance methodology, the engagement sequence of a short-range, 
air-defense, missile system is described. The Roland system 
was selected because the general operational aspects of the 
system are available at the unclassified level [1] . 

The entry of one or more aerial targets into the range of 
the search radar is indicated to the Roland vehicle commander 
by an audible tone . At the same time , a synthetic display of 
the targets appears on a screen to give the commander all the 
information needed to select the most threatening target. The 
screen images are different for friendly and enemy targets. 
Also, the entry of the target into the missile envelope, 
utilizing target advanced-range computations, is indicated by 
a change in the display . With the search antenna raised and 
the search radar activated, target acquisition is possible even 
when the vehicle is in motion. 

There are three modes of identification, friend or foe 
(IFF) interrogation: automatic, manual, and automatic within 

a given range . 

When the commander has recognized a target as hostile 
and decided to engage it, he places a cursor over the screen 
image. This automatically brings the turret to bear and 
tracking can commence in either the "radar" or "optical" modes. 
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In the "radar" mode, the tracking radar automatically 
accepts target designation from the search radar, searches 
for, locks onto, and tracks the target. 

In the "optical" mode, the aimer searches for the target 
in elevation with an optical sight. To aid him an electronic 
instrument displays the maximum theoretical elevation for 
the search. When the aimer has acquired the target in his 

J cross-hairs, he keeps the target in his sight by manipulating 

B 

a control stick. This control keeps the target properly 
positioned by moving the turret in azimuth and swivelling a 
mirror in elevation. 

As soon as the commander confirms that the target is 
within missile range, he initiates the firing sequence in 
the "radar" mode, or authorizes "optical" mode firing through 
a command displayed in the aimer's sight. The aimer, then, 
can initiate the firing sequence. 

The missile is guided by a command to line-of-sight 
technique. This means that the target is tracked optically 
or by radar and the deviation of the missile from this line 
of sight is determined and corrected by a guidance command. 

; The commander may switch from "radar" to "optical" and back 

[ again, as desired, even after the missile has been launched. 

I 

I Target tracking and determination of the missile's deviation 

1 

y from the line of sight are different for each mode. In the 

"radar" mode, the guidance radar has two receiving channels. 

One is used for target tracking and the other is used to 
locate the missile in the radar lobe through reception of the 
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missile's radio frequency beacons. By comparing these angles, 
an error between the missile and the target line of sight can 
be determined. In the "optical" mode, a biaxially-stabilized 
mirror is manually controlled to keep the target vertically 
in the aimer ' s sight and the turret is rotated to the azimuth 
j of the target line of sight. An infrared goniometer is mounted 
i to provide missile angle from the tracker by following flares 
ti mounted on the rear of the missile. Then, a deviation of the 
missile angle from the target line of sight can be determined. 

Two groups of signals are introduced into the command 
computer: the velocity of the line of sight in azimuth and 

elevation, and the deviation of the missile from the line of 
sight in azimuth and elevation. Based upon data from the line- 
of-sight movement and the angular deviations of the missile, 
the necessary guidance signals are calculated. 

The guidance signals are relayed to the missile by a 
transmitter with highly directional characteristics. The 
command- transmitting antenna is slaved to the missile angle 
in both azimuth and elevation. It, therefore, is trained 
on the missile continuously. 

The side forces required for missile course corrections 
are produced through deflection of the exhaust jet of the 
’ sustainer motor by spoilers at the rear of the missile (thrust- 

I 

ij vector control) . 

When the missile reaches the point of impact with the 
target, the warhead is detonated by either percussion, contact 
fuse or the radio- frequency, proximity fuse. The warhead 
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consists of a radial-effect, multiple- fragmentation charge. 

Figure III-l presents an operational schematic of the 
basic Roland missile system operation. 




Figure III-l Roland Missile System Operational Schematic 



The computer simulations contained herein are generic in 
nature within the command to line-of-sight guided-missile type 
and have only reasonable estimates of missile capabilities 
introduced. This ensures unclassified results. At the same 
time, the simulations are of sufficient complexity to 
properly weigh the relative merits of the guidance variations 
discussed . 
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IV. COMMAND TO LINE-OF-SIGHT GUIDANCE GEOMETRY 



There are many tactical, inherently short-ranged, surface- 
to-air missiles which use this guidance technique for their 
operation. Its objective, basically, is to keep the error 
term (e) zero so that the missile flies directly up the line 
of sight. 

Since only angles to the target and to the missile are 
required for guidance, attempts to preclude target range 
through denial jamming prove useless against this type of 
guidance system. The only necessary information is, in general, 
to know whether the target is within the missile envelope. 
Deceptive down- link jamming, however, appears to be feasible. 

One possible method is presented in Section VI . 

The geometry depicted in Figure IV-1 summarizes the 
"basic" command to line-of-sight geometric situation. 




Figure IV-1 "Basic" Command to Line-of-Sight Geometry 
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A variant of this "basic" conunand to line-of -sight guidance 
approach, which shows promise for increasing system effectiveness, 
is to incorporate "lead-angle" into the basic guidance model. 

This technique would increase the effective range of the missile 
against a non-maneuvering, crossing target; however, it is de- 
graded by a highly-maneuvering target. Since range information 
is required to calculate the impact point, additional electronic 
countermeasures could be employed by the target to deny this 
range information thus negating the "lead-angle" command to line- 
of-sight missile guidance scheme. 

The previous "basic" geometry has been retained but it has 
been modified to utilize a synthetic line of sight upon which 
the missile flies. This synthetic line of sight is the direction 
of the impact point from the ground tracking unit. Therefore, 
the missile is commanded to fly toward the impact point rather 
than directly toward the target. This method, then, reduces 
the missile's flight time until impact, thus, improving the 
missile system's effective engagement range. Figure IV-2 shows 




Figure IV-2 "Lead-Angle" Command to Line-of-Sight Geometry 
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V. ON-OFF, THRUST- VECTOR, MISSILE CONTROL 



The guidance models used in the computer simulations utilize 
on-off control. This type of control also is known as relay 
control or "bang-bang" control because the system's actuator is 
not moved in proportion to the error signal but, rather, is 
positioned either totally in the positive direction or totally 
in the negative direction. This causes the missile's thrust- 
vector control to apply missile-maneuvering thrust fully in one 
direction or the other. At first glance this type of system 
might appear to be very crude; however, it does have major advan- 
tages. First, the system does not require complex proportional 
amplification; the unit is smaller than its proportional control 
counterpart; it is, generally, the time-optimal control system; 
it weighs less than other systems; and, finally, it is relatively 




Figure V-1 Missile Cross-Range Error with Zero Initial 
Values for Cross-Range Error and Error Rate 
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Based upon this type of missile control system, initial 
computer simulations, not presented here in detail, were designed 
They provided an elementary computer program structure which was 
modified incrementally to generate improved computer models. In 
these first programs , the error values were introduced directly 
into the missile's thrust-vector control without error rate. 

This produced increased missile cross-range error with flight 
time as shown in Figure V-1. This effect has been eliminated 
in the subsequent computer guidance models by proper acceleration 
limiting techniques. Additionally, if initial erroneous values 
of cross-range error and cross-range error rate are introduced 




Cross-Range Error (m) 

Figure V-2 Phase Plane When Control Function By-Passed 

With Initial Values Introduced for Cross -Range 
Error and Error Rate 
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just after missile launch, this "basic" guidance system does 
not dampen their effects. This increases the probability of 
having a greater miss distance. With an initial cross-range 
error of ten meters and with an initial cross-range error rate 
of fifty meters per second applied, the undamped effect, plus 
increased cross-range error with time, combine to cause the 
' missile to miss the target. Results are shown in Figures V-2 
i through V-4. Although not readily apparent in Figure V-4 , the 




Time (sec) 



Figure V-3 Missile & Target Positional Coordinates When 

Control Function By-Passed With Initial Values 
Introduced for Cross-Range Error and Error Rate 
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X-Position (ra) 



Figure V-4 Missile & Target Flight Paths When Control 
Function By-Passed With Initial Values 
Introduced for Cross -Range Error and Error 
Rate 



Cross-Range 

Error 

(m) 




Figure V-5 Missile Cross-Range Error When Control Function 
By-Passed With Initial Values Introduced for 
Cross-Range Error and Error Rate 
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missile moves about the line of sight with a cross-range error 
varying as shown in Figure V-5. The cross-range error does 
increase by a small amount with time as shown in Figure V-1. 

The models discussed later are somewhat unrealistic because 
they assume that the required force obtained from the missile's 
thrust-vector control can be applied instantaneously at the 
proper time and reversed in a similar fashion to decelerate the 
missile. This, however, does provide optimum missile performance 
for comparison with any real system. 

In all the simulations the missile is acceleration-limited 
with instantaneous application of constant side-force assumed. 

One could allow for a time delay in switching the missile ’ s 
thrust-vector control from the positive to the negative direction. 
However, this would raise the order of the analysis. The computer 
program, though, can be modified to handle systems to the 
thirtieth order. 

Optimum acceleration switching requires that a switching- 
boundary be established that zeros the error in the minimum time. 
That boundary was determined to be parabolic as drawn in the 
phase-plane representation of Figure V-6 and as supported 



CRE 
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by the calculations which follow later. When the acceleration 
or deceleration trajectories attempt to cross this optimum 
curve (switching boundary) the system switches to that optimum 
path and the error is driven to zero along it. 

! In order to verify this result, a phase plane was produced 

by introducing, as initial conditions, large cross-range error 
I and cross-range error rates into a computer-simulated guidance 
ii model containing acceleration-limiting methods within the tracker 
control function. The resulting phase plane. Figure V-7 , sub- 
stantiates that the tracker control function switches properly 
to the theoretically-optimum, parabolic curve for error reduction 
expected by calculation. This curve is represented in Figure 
V-6. Also, the target was intercepted by the missile as sub- 
stantiated in Figures V-8 and V-9. Figure V-10 shows the 



Cross -Range 
Error Rate 
(m/sec) 



!l 







Cross-Range Error (m) 



Figure V-7 Phase Plane When Control Function Utilized 
With Initial Errors Introduced for Cross- 
Range Error and Error Rate 
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Coordinate 

(m) 




Time (sec) 



Figure V-8 Missile & Target Positional Coordinates When 
Control Function Utilized With Initial Errors 
Introduced for Cross-Range Error and Error Rate 



Y-Position 

(m) 




Figure V-9 Missile & Target Flight Paths When Control 
Function Utilized With Initial Errors 
Introduced for Cross-Range Error and Error 
Rate 
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Time (sec) 

Figure V-10 Missile Cross-Range Error When Control 
Function Utilized With Initial Errors 
Introduced for Cross-Range Error and 
Error Rate 



missile's cross-range error with time of flight. This error is 
not readily apparent in Figure V-9 because of the relatively 
large distance scale utilized in that figure. As can be seen 
this cross-range error is reduced quickly. 

The theoretical calculations supporting the parabolic 
nature of the switching boundary represented by Figure V-6 
follow: 

m = missile mass 

F = lateral force applied to missile 

CRE = cross-range error 

CRE = cross-range, error rate 

• • 

CRE = cross-range error acceleration 

G = magnitude of cross-range missile acceleration 
U = cross-range missile acceleration 
C = constant 

Since the mass of the missile and the magnitude of 
the output force from the missile thrust-vector 
control are assumed constant, the magnitude of the 
acceleration is constant, thus: 

• • 

m CRE = |f| 

• • 

I CRE I = |F|/m = G 

I - CRE I = G 
dt 

J* I d CRE I = Jg dt ^1 CRE | = Gt + c 

Let the initial conditions for time and cross-range 
error be zero. Therefore, the constant, "C", 
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is zero also. 



CRE = Gt , CRE - 0 
t = CRE/G ^ = CRE^/G^ 

I CRE = Gt J*d(CRE) = jGt dt 
CRE = Gt^/2 + |CRE(0) | , = CRE^/G^ 

CRE = CRE^/2G + |CRE(0) | , CRE - 0 

Since the direction of the missile's thrust-vector 
control force can be either in the positive or 
negative direction, the sign of the cross-range 
error (CRE) must be preserved. The initial missile 
cross-range error can be in either direction; 
consequently, the equation becomes: 

CRE = CRE |CRE|/2G + CRE(O) 

This equation is parabolic as depicted in Figure 

V-6. 

Also, since the magnitude of the missile cross- 
range acceleration, "G" , is constant, the correspon- 
ding acceleration vector, "U" , is: 

U = ±G j U = (-G) SIGN (error function) 

, error function = EX = CRE | CRE | /2G + CRE 

This drives the missile's control surfaces in the proper 
direction, a direction which is opposite to the sign 
of the error function, thus moving the missile back 
toward the line of sight with a constant magnitude of 
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acceleration, ”G". 



A simplified block diagram representing the command to 
line-of-sight guidance system is presented in Figure V-11. 




Figure V-11 Elementary Command to Line-of -Sight Block 
Diagram 
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VI. 'BASIC COMMAND TO LINE-OF-SIGHT SIMULATION 



A. PROGRAM DESCRIPTION 

Utilizing the concepts previously discussed, this program 
simulates theoretically optimum surface-to-air missile maneuver 
against a constant-velocity target in two dimensions. The 

li engagement was designed to occur within the first quadrant of a 

'i 

i standard Cartesian-coordinate reference system with the ground 

' tracker, missile lauching unit, located at the origin. The 
missile was assumed to be "captured" by the guidance system 
approximately one second after missile launch or at a distance 
of about 500 meters from the ground tracker unit. At the time 
of missile "capture" an initial cross-range error of ten meters 
and a cross-range error rate of 50 meters per second were intro- 
duced to demonstrate the ability of the guidance system to neu- 
tralize initial errors. It should be noted that these values 
of error were selected entirely arbitrarily and that they could 
be any reasonable values assumed to exist at the time of missile 
"capture" . 

The target was "flown" across the first quadrant from a 
position 5000 meters from the tracker on the X-axis (5000,0) 
toward a position 5000 meters from the tracker on the Y-axis 
(0,5000). 

I In order to use a command to line-of -sight guidance scheme 

effectively, the missile must possess a speed advantage 
over the target. Therefore, in this model the missile velocity 
was chosen to, be 500 meters per second versus a 250 meter-per- 
second target. At sea level these speeds equate to approximately 
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Mach 1.5 and Mach 0.75 respectively. The missile, also, after 
being launched directly toward the target, was maneuvered using 
on-off, thrust-vector control which applied a lateral acceleration 
of 60 meters per second squared to the missile. This is just 
greater than six times the earth's gravitational acceleration 
at sea level. 

In general, keeping in mind the scenario above, the program 
operates as follows : 

The equations for the target ' s movement were 
written and appropriately placed into the program 
so that the target and missile angular directions 
could be compared. The difference between these 
two angles creates an error angle. This error 
angle is geometrically converted by the tracker 
computer into the cross-range error, the distance 
of the missile from its intended flight path - 
the line of sight to the target in this case. 

That error distance is reduced to zero by 
applying a lateral, thrust-vector-controlled 
acceleration to the missile in the correct 
direction as determined and commanded by the 
ground tracker unit. The optimum parabolic 
switching boundary, discussed in Section V, was 
used in the tracker computations . In order to 
determine this cross-range error , the tracker 
estimates the missile's range by the elapsed time 
of flight and the missile's velocity profile. 

Also, to simplify the program somewhat, the 
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missile's velocity was assiamed to be constant 
in a radial direction from the tracker. When the 
cross-range error is reduced to an acceptable dis- 
tance, the amount of thrust-vector control is 
limited to prevent rapid changes between full 
thrust-vector deflections in both directions. If 
the error, again, becomes unacceptably large, the 
full deflection positions are utilized. Missile 
position is determined by state variable analysis 
using the INTEG2 subroutine for integration of the 
lateral thrust-vector-controlled missile acceleration. 
The use of this subroutine is explained in Appendix 
A. After two integrations, the missile's new cross- 
range error is determined and, then, this distance is 
used with estimated missile range to provide a revised 
missile angular position. This process is repeated 
in increments of five thousandth of a second from 
the time of missile "capture" by the tracker, a 
time of approximately one second, until maximum 
missile range is attained after ten seconds of 
flight. Thus, the missile is driven continually 
to a position upon the line of sight to the target 
The program variables are "frozen" when the missile 
reaches the closest point of approach to the target 
if the value of constant C (1) is one. If this 
value is zero, the program computes until the time 
of missile flight reaches ten seconds. 
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The program was written in the FORTRAN IV language. A 
detailed, annotated copy of the program, including a table of 
program variable definitions, is provided following the 
narrative portion of this paper. 




Figure VI-1 'Basic' Computer Simulation Block Diagram 
B. SYSTEM BLOCK DIAGRAM 

The command to line-of-sight computer simulation block 
diagram. Figure VI-1, shows how the elementary block diagram 
of this guidance technique. Figure V-11, was expanded upon in 
order to develop a computer model for the system. As will be 
seen, while appearing to be quite simple in concept, the 
simulation becomes increasingly more complicated as its sophis- 
tication is enhanced. 

C- SIMULATION RESULTS 

The plotted output of the phase plane. Figure VI-2, shows 
how the cross-range error and cross-range error rate are driven 
to a near zero condition along the optimum, parabolic, switching 
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boundary. Once the cross-range error is less than one meter, 
the magnitude of the thrust vector control is reduced by a 
sixth to minimize the possibility of over-controlling the 
missile. 



Cross -Range 
Error Rate 
(m/sec) 




Cross-Range Error (m) 

. Figure VI-2 'Basic' Simulation Phase Plane 



Figure VI-3 presents the "X" and "Y" positions for the 
missile and the target versus time. As can be seen the inter- 
cept occurs in just under eight seconds after missile launch. 

Figure VI-4 indicates the relative positions of the 
missile and target with the terminal points for each flight 
profile stopped at the closest point of approach. An intercept 
did occur. 

In order to ensure that the missile's cross-range error 
from the line of sight was damped. Figure VI-5 was produced to 
show this cross-range error versus time of flight. The initial 
errors, applied at the time of missile "capture" by the 
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Figure VI-3 'Basic' Simulation Missile & Target 
Positional Coordinates 
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Figure VI-4 ' Basic' Simulation Missile & Target 
Flight Paths 
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Figure VI-5 'Basic' Simulation Missile Cross-Range Error 
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guidance system (corresponding to a time of approximately one 
second in the figure) , were reduced sufficiently in less than 
two seconds . 

D. EFFECTS OF DECEPTIVE JAMMING 

A readily apparent deceptive jamming technique, which might 

j be effective against this class of guidance system, would be to 

j employ, aboard the target aircraft, an emitter which imitates 

' the signal that is used by the tracker to locate the missile's 

1 

angle from the reference plane. This could be, for example, 

1 ' 

depending upon how the missile is tracked in angle, an imitation 
of the missile's transponder radio-frequency signal or a dup- 
lication of the missile's infrared emissions. In any case the 
deceptive jamming signal from the target aircraft should be of 
sufficient strength to ensure that the missile tracking system 
within the ground- located tracker considers the target to be 
the missile. There are, most assuredly, counter-countermeasures 
which could be applied to negate this deceptive technique; how- 
ever, this, currently, would be one method of generating a zero 
error signal between the missile and target angles. Additionally, 
this deceptive countermeasure provides a method to demonstrate 
how the computer simulation drives the missile with inputs of 
erroneous guidance information. 

^ In order to implement this deceptive jamming technique, it 

I 

j was assiimed that the target aircraft jamming equipment became 
effective when the missile was approximately half the distance 
to the target. This would equate to about four seconds of 
missile flight time. Therefore, since this jamming 
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countermeasure system simply zeros the error angle, only one 
additional equation was necessary to modify the previously 
discussed "basic" command to line-of-sight guidance model. The 
equation, IF (T.GE . 4 . 0) E=0 . , was placed into the program following 
the computational step for the calculation of the angular error, 
E=SIGM-SIGT. 

The resulting missile trajectory and other parameters are 

[;| readily apparent in Figure VI-6 through VI-8. As shown in 

p 

; Figure VI-6, the error was reduced to zero by the time that 

i' 

I' jamming commenced. After deceptive jamming was initiated, the 
cross-range error rapidly increased causing a miss of over 300 
meters. Needless to say, because of the short times of flight 
involved, the countermeasures must be employed in a timely 
manner. Otherwise, the missile could approach, easily, within 
lethal range of the target aircraft. 
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Figure VI-6 ' Basic' Simulation Phase Plane When 
System Deceptively Jammed 
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Figure VI-7 'Basic' Simulation Missile & Target 
Flight Paths When System Deceptively 
Jammed 




Figure VI-8 'Basic' Simulation Missile Cross-Range Error 
When System Deceptively Jammed 
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VII. 



' LEi^D- ANGLE' COMMAND TO LINE-OF-SIGHT SIMULATION 



A. PROGRAM DESCRIPTION 

This 'lead-angle' simulation uses the same target profile as 
that of the 'basic' command to line-of-sight simulation. This 

I allows a comparison of the two simulations to a certain degree. 

!1 

One parameter, which can not be compared accurately, is the time 
to impact. This occurs because the radial missile velocity from 
the tracker was assumed to be a constant value - in this case 
500 meters per second. Therefore, intercepts will take place 
at about the same time regardless of missile maneuvers. A 
general comparison can be made, however, by mentally comparing 
the plotted missile flight profiles, the longer profile taking 
the most time in an actual engagement. Also, in order to avoid 
redundancy, discussion will be limited to modifications of the 
previously described computer model. Only a few changes were 
required to produce a 'lead-angle' variant for the command to 
line-of-sight guidance system. 

In this guidance scheme, the special positions of both the 
target and the missile must be accurately and continuously 
known in order to compute an advanced impact position. This is 
accomplished by calculating the distance from the missile to 
the target and their closing velocity. From this data, a time 
j remaining until impact can be estimated . This remaining time 
i of flight, then, is combined with the target's current velocity 
I vector and position to ascertain the expected impact point. A 
synthetic line of sight toward this location is generated by 
the ground tracker unit computer and the missile is commanded 
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to "fly" up this flight path instead of along the line of sight 
directly toward the target. The guidance scheme used to main- 
tain the missile on course, the synthetic line of sight, is, 
otherwise, identical to that of the "basic" simulation previously 
described. Once the missile passes the closest point of 
approach to the target, this simulation ceases to be valid. This 
is quite all right because, once the missile has passed the 
target, no maneuver would be sufficient to effectively engage 
the aircraft. 

B. SYSTEM BLOCK DIAGRAM 

The "lead-angle" command to line of sight computer simulation 
block diagram is contained in Figure VII-1. It is essentially 
a modification of the "basic" simulation block diagram. Figure 
VI-1. As described previously, target and missile positions 
plus their velocities must be known for the guidance scheme to 
be viable. 
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c. SIMULATION RESULTS 

This simulation supports the idea that the use of a "lead- 
angle” variant to the "basic" command to line-of-sight guidance 
system enhances the effectiveness of the elementary guidance 
technique. As can be seen in Figures VII-2 through VII-5, the 
missile was "captured" by the guidance system at a time of 
approximately one second from launch. The initial cross-range 
and cross-range-rate errors were quickly reduced and the missile 
"flew" out the synthetic line-of-sight to successfully intercept 
the target aircraft. Of course, if the target employs electronic 
countermeasures such as denial jamming of the target tracking 
equipment, this "lead-angle" method would prove to be useless 
because accurate positioning information must be maintained to 
compute the estimated impact point. Also, if the target 
maneuvers sufficiently to cause the impact point to vary at 
high rates, the synthetic line of sight may shift at too high 
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Figure VII-2 'Lead-Angle' Simulation Phase Plane 
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a rate for the missile to follow. Thus, the missile may not 
be able to physically engage the impact point but might be 
able to still effectively attack the target by the "basic" 
guidance scheme. Consequently, in order to capitalize upon 
the advantage of the "lead-angle" variant (an increased effec- 
tive missile envelope against a relatively slow maneuvering 
target which is not employing electronic countermeasures) , 
it seemed that a combined guidance technique should be con- 
sidered. Therefore, the "consolidated" simulation discussed 
in the next section was contrived. 
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Figure VII-5 'Lead-Angle' Simulation Missile Cross-Range 
Error 
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VIII. 



' CONSOLIDATED ' COMMAND TO LINE-OF-SIGHT SIMULATION 



A. PROGRAM DESCRIPTION 

When the two previous simulations were combined, a major 
problem surfaced. That problem was basing the missile's new 
position upon two different references , the line of sight to 
the target in the case of the "basic" model and the synthetic 
line of sight to the impact point in the case of the "lead-angle" 
variant. This incompatibility caused the missile's position 
to jump radically as the "consolidated" simulation attempted 
to shift alternately between the two reference lines as the 
jamming situation varied. Therefore, a single reference had 
to be implemented for continuous use regardless of whether 
jamming was or was not encountered. The X-axis of the Cartesian- 
coordinate system was selected since various angular measurements 
were based upon it already. Also, this has proved to provide 
a more accurate model of the missile's optimum flight profile 
because the reference does not move during the time of missile 
flight. If a strictly "basic" simulation should be desired, 
jamming is considered to be effective during the entire engage- 
ment. Conversely, the absence of jamming during the missile's 
flight produces a total "lead- angle" intercept of the target 
aircraft. Of course, different duty cycles for jamming can be 
implemented by the target during a single computer run to test 
the ability of the missile system to transition between the two 
attack modes. A limited number of different jamming combinations 
were tested and are discussed later. 
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In order to implement jamming, the jamming time increments 
are read into the simulation through five time intervals 
specified on the first five data cards- These cards preceed 
those required by the INTEG2 subroutine as explained in 
Appendix A- The time each jamming interval begins is punched 
into columns 1-10 of each card in F-format. The times that 
jamming ceases are entered in columns 11-20. Should a no- 
jamming condition be desired, then, all five cards are left 
blank but they must not be ommitted. Total jamming can be 
obtained by making one card cover the entire time of flight. 
For example, the first card could have colxomns 1-10 blank and 
the value 11.0 entered into columns 11-20- This produces 
jamming from time zero until a time of eleven seconds, which 
is well after the end of the engagement. The remaining four 
cards are blank and must be inserted. 

This simulation "drives" the missile in the "Y" direction 

with sufficient "Y" velocity to reduce the cross-range error 

to zero. The missile's X-direction velocity component is used 

to maintain a constant, outward, radial velocity, "V ", toward 

m 

the target. In this simulation, as in its predecessors, that 
radial velocity, , is assumed to be 500 meters per second. 

Therefore, the velocity of the missile in the "X" direction is 
related to that obtained from the forcing function in the "Y" 
direction. This is shown in Figure VIII-1. 




Figure VIII-1 'Consolidated' Simulation Missile Velocity 
Vectors 



48 



The missile's Y-direction velocity is obtained from the state 
variable analysis through a single integration of the forcing 
function, "U" , which is the missile's Y-direction acceleration. 
This acceleration is the "Y" component of the cross-range error 
acceleration achieved by thrust-vector control of the missile as 
shown in Figure VIII-2. 




Figure VIII-2 'Consolidated' Simulation Missile 
Acceleration Vectors 

The missile's Y-coordinate position, then, is the second inte- 
gration of the forcing function. From these results, the X- 
direction components were determined since the velocities must 
yield a constant, radial velocity of 500 meters per second. In 
order to ensure a correct integration process, the missile's 
"Y" velocity was not limited to 500 meters per second; however, 
it never exceeded that value in any of the computer runs tested. 
It should be noted that the program will zero the acceleration 
and the velocity of the missile in the "X" direction if excessive 
missile "Y" velocity is experienced. Also, the results, then, 
will become invalid; consequently, the velocity of the missile. 



49 



i 



i 



1 



"V ” , should be increased or the target ' s speed reduced to 
m 

obtain a theoretically-possible engagement scenario. 

The targets, in the previous simulations, "flew" a constant- 
velocity, linear course. In this simulation the target was 
maneuvered sinusoidally in the "Y" direction for the purpose 
of observing the effects of a moving impact point. The maximum 
velocity of the target in the "Y" direction is 64 meters per 
second with a half period of four seconds. The target's "X" 
velocity was a constant, negative 210 meters per second. This 
did produce some interesting results. First, this type of tar- 
get maneuver caused the impact position to move at even higher 
rates than the target moved. Thus, in some cases, the missile 
may not be able to engage, adequately, the impact point which 
is moving at a very high velocity. It, however, may be able to 
attack the target itself by utilizing the "basic" command to 
line-of-sight approach. To determine which situation exists 
necessitated the use of the impact point's velocity to compute 
the synthetic line of sight, the angle upon which the missile 
"flies" in a non- jamming environment. If this impact-point 
velocity exceeded eight-tenths of the missile's velocity, then, 
the missile was commanded to "fly" directly toward the target's 
present position by the "basic" guidance method. This was in 
lieu of utilizing the "lead-angle" approach with only an un- 
stable impact position and, thus, an erratic or too quickly 
moving synthetic line of sight. 

In order to provide a longer time of flight, the target's 
initial position was moved to coordinates (6000,1000) and the 
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target was "flown" inbound at 210 meters per second while 
sinusoidally maneuvering about the 1000 meter Y-ordinate. Also, 
in regard to the short time of flight involved, excessively 
large initial errors in either distance or rate should be 
avoided because of the time required to neutralize these errors . 
For example, since the missile's maximum acceleration is only 
150 meters per second squared, initial erroneous conditions of 
250 meters and 250 meters per second for the missile's Y-direction 
components require many seconds to correct. However, with 
initial conditions of 150 meters and 150 meters per second, the 
simulation proved the missile to be effective against this 
target over a full scope of jamming duty cycle variations - 
from none to continuous denial jamming. 

The initial angle to the missile at the time of "capture" 
by the guidance system, a time of one second or 500 meters from 
the launcher-tracking unit, is determined by the initial con- 
ditions, X(l) and X(2). As can be recalled, X(l) is the 
missile's "Y" position at the time that the weapon is 500 meters 
radially from the launcher. Therefore, the angle is the inverse 
sine of X(l) divided by 500 meters. In order to attain this "Y" 
position in one second, the missile's "Y" velocity, then, was 
approximately the "Y" distance per second. Consequently, in 
this simulation, initial conditions of 100 meters for X(l) and 
100 meters per second for X{2) equate to an angle of 11.54 
degrees above reference. 

Also, the magnitude of the missile's maximum lateral 
acceleration, "G", was increased to 150 meters per second 
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squared. This provided a sufficient forcing function to allow 
engagement of the target. This was possible even with increased 
initial errors at the time of missile "capture" by the guidance 
system; with a maneuvering target accelerating at just over 50 
meters per second squared, at maximum, in the "Y" direction; 
and with various jamming duty cycles from zero to one. As in 
the previous simulations, the thrust-vector control of the 
missile was reduced by one-sixth; however, this was done only 
when the sum of the absolute values of cross-range error and 
error rate was less than or equal to four vice the value of one, 
which was used previously. The error tolerance was increased to 
provide a larger "idling" area for the missile's thrust-vector 
control because of the increased forcing function utilized. 

The actual computer program for the simulation, with 
comments included, follows the appendices. 

B. SYSTEM BLOCK DIAGRAM 

The block diagram of the "consolidated" simulation is con- 
tained in Figure VIII-3 . Although it appears formidable at 
first glance, the diagram merely combines the concepts of the 
two previous simulations with slight modifications to account 
for the use of a new reference from which to compute the 
missile's location. If the tracking unit ascertains that it is 
contending with either denial jamming or a highly-maneuvering 
target, then, only the upper portion of the tracker diagram is 
used for guidance.. Conversely, if both a non- jamming condition 
and a relatively slowly -maneuvering target situation are 
experienced, the lower portion of the tracker diagram produces 
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guidance for the missile. 

The simulation does produce necessary calculations in all 
areas of the tracker diagram continually throughout the engage- 
ment, but only generates commands to the missile from the 
appropriate section of the model. 

After the program reads the data cards , the times of 
jamming are known and the appropriate control function within 
the tracker's computer is implemented to control the missile. 

The magnitude of the impact point's velocity is considered, also, 
to ascertain the stability of the synthetic line of sight and, 
thus, to determine whether target maneuvers make the "lead- 
angle" approach of guidance inappropriate. 




Figure VIII-3 'Consolidated' Computer Simulation Block 
Diagram 
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C. SIMULATION RESULTS 



Three engagements were computed under different jamming 
conditions. They were continuous jamming, intermittent jamming 
having an approximate duty cycle of one-half, and no jamming. 
During intermittent jamming, the ground- located tracker unit 
was denied range information within the following time increments 
after missile launch: 1. 0-3.0, 4. 0-5.0, 6. 6-7.0, 7. 4-7. 8, and 

8. 2-8. 6 seconds. 

Utilizing the scenario which was previously discussed, all 
three attacks resulted in successful intercepts of the target. 

As can be observed in Figure VIII-4 , the best flight profile 
of the three simulated runs occurred when total jamming was 
experienced. This jamming caused the missile to "fly" 
directly toward the target. 




Figure VIII-4 'Consolidated' Simulation Missile & Target 
Flight Paths When Continuously Jammed 
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Figure VIII-5 'Consolidated' Simulation Missile & Target 
Flight Paths When Not Jammed 




Figure VIII-6 'Consolidated' Simulation Missile & Target 
Flight Paths When Intermittently Jammed 
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When no jamming was encountered, the flight profile became 
more exaggerated because the target was maneuvering , thus 
causing the estimated impact point to move at an even higher 
rate. Recall that, in a non- jamming environment, the missile 
will "fly" for the impact point vice the target. Figure VIII-5 
shows this profile. Finally, with intermittent jamming, the 
missile's flight profile. Figure VIII-6, became more erratic 
but did not reach the more extreme deviations experienced in the 
un jammed case. Therefore, it appears that the best method would 
be to use the "basic" command to line-of-sight guidance technique 
against a maneuvering target regardless of the jamming condition 
in effect. Also, as can be observed in Figure VIII-6, the flight 
profile during intermittent jamming, the amount of erratic 
variation in the missile's movement decreased considerably as it 
closed the target. This occurred because the difference between 
the synthetic line of sight and the line of sight becomes 
smaller with time until they are equal in value at the time of 
impact. Consequently, the "lead-angle" approach should be 
considered for utilization against a non-maneuvering target 
when adequate range information is available, i.e., when denial 
jamming is not encountered. 

In order to reduce the amount of movement in the missile's 
thrust-vector controllers and, thus, to ensure that they do 
not "beat themselves to death" mechanically while trying to 
reduce a trivial cross-range error, a reduced amount of 
deflection was applied when the cross-range error was relatively 
small. This, then, was a trade-off of allowing some increase 
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in cross-range error for improved thrust-vector control 
mechanical reliability. Figure VIII-7 shows how the cross- 
range error was reduced without using a smaller amount of thrust- 
vector control deflection. Keep in mind that the control sur- 
faces in the missile were shifting rapidly from a total positive 
to a total negative direction. In Figure VIII-8, a cross- 
range error of four meters was accepted as a reasonable maximum 
miss distance and the amount of control was decreased when the 
error was within that value. This assured improved mechanical 
reliability of the system. As can be seen in Figure VIII-8, 
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at the time of impact the cross-range error was in the process 
of being reduced from a value of about four meters. The phase 
planes, Figures VIII-9 and VIII-10, also, show the effects of 
this trade-off. Figure VIII-9 was the result obtained when 
using no reduction in the amount of thrust-vector control; 
while Figure VIII-10 shows the effect of allowing a maximum 
four-meter error with a reduced amount of lateral thrust. 



IX. CONCLUSIONS 



These simulations support the premise that consideration 
should be given to utilizing a "lead-angle" variant in con- 
junction with the "basic" command to line-of-sight guidance 

\ 

I technique. This would increase the missile system's effective 
engagement envelope against a non-maneuvering target that 
employs low- duty- cycle denial jamming. If the target utilizes 
higher duty cycle jamming or begins to make evasive maneuvers, 
then, the system can revert, when necessary, to the "basic" 
guidance mode even after the missile has been launched. 

This new approach appears to be genetically applicable 
to missiles of the command to line-of-sight guidance class. 

If a range- tracking capability already is contained within 
the system, it should be possible to implement this improved 
guidancfe technique through modifications to the existing 
guidance computer. 
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APPENDIX A 



DESCRIPTION AND USE OF SUBROUTINE INTEG2 



INTEG2 provides a simple means for a Runga-Kutta solution 
of ordinary differential equations. Input and output routines 
(including graphical point plots on the versatec plotter) are 
built in. 

The subroutine has been compiled and added to the FORTRAN 
library, SYSl.MPSLIB. The program is flexible and self con- 
tained. It will seldom be desirable to build the subroutine 
into a more complex program. 

It is assumed that the equations to be solved have been 
reduced to 'N' simultaneous first order differential equations 
of the form: 

DX(I)/DT = F (X(l) ,X(2) , . . . ,X(N) ,T) , I = 1,2, ...,N 
where the function 'F' may be nonlinear. The user must trans- 
late these 'N' equations into FORTRAN as in the following second- 
order example: 

XDOT(l) = X(2) 

XD0T(2) = C(l) *X(1)-0.2*X(2)+C(2) *SIN(T*T) 

C(l) and C(2) are constants which are entered on data cards. 

To solve these equations the OS/360 deck must be of the 
following form: 

// STANDARD JOB CARD (GREEN) 

// EXEC FORTCLGV, REGION. G0=150K 

//FORT. SYS IN DD * 

DIMENSION X(30) ,XDOT (30) ,C(15) 
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C(10) = 1. 

1 CALL INTEG2 (T,X,XDOT,C) 

• • • 

USERS FORTRAN EQUATIONS 
• • • 

GO TO 1 
END 
/* 

//GO. SYS IN DD * 

DATA CARDS AS DESCRIBED BELOW 
/* (ORANGE) 

The number of XDOT equations is 'N', which must not exceed 
30. These equations are supplied by the user, who defines 
each 'XDOT' in terms of the dependent variables X(l) through 
X(N) and the independent variable 'T'. In writing these 
equations the user may introduce at his convenience: 
any unsubscripted variables; 

the subscripted variables X(I) , where the variable 
'I' must be less than or equal to 30; 
the constants C(I) to be entered as data, where the 
variable 'I' must be less than or equal to eight; 
any normal FORTRAN technique or function; and 
routines from any source library or user- supplied 
subroutines . 

Note 1: The use of 'auxiliary' X(I), when the 

variable 'I' is greater than the number of equations 
to be integrated, does not alter the value of 'N' , 
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the order of the equations. 

Note 2: Loops, either with or without a 'DO' state- 

ment, are best avoided. 

Note 3 : ' IF ' statements , provided that they do not 

create a loop, can be used to transfer control 

i’ 

I within the user's equations. For example, the 

, statement, IF (T. GT . 10 . ) C (3) =0 . , would cause C(3) 

I 

to take the value zero for all 'T' greater than ten. 

« 

The constants C(l) through C(8) may be used as described in 
the above examples, and are read in from a data card. C(10) 
must never be used, except as indicated in the standard OS/3 60 
deck above. C(ll) and C(12) control the output. For example, 
if the statements, C(ll) = 10.0 and C(12) = 2.0, are added to 
the users FORTRAN equations, every tenth integration step will 
be printed out, and every second step will be plotted. (If 

S 

not set by the user, default values of 20 and five apply.) 

C(13) can be similarly used to modify the step size of the 
numerical integration which is more usually defined on a data 
card. 

Only the independent variable 'T' and the variables X(l) 
through X(30) can be output. Therefore, to output a quantity 
which is not one of the 'N' dependent variables. The user 
must introduce an 'auxiliary' variable X(I), where the variable 
I 'I' must be less than 31 but greater than 'N' - for example, 

: by introducing the FORTRAN equation, X(27) = X(3)*X(3), the 

square of X(3) can be output by the program as X(27). Normally, 
a line of printout is generated after every 20 integration 
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steps, and one graph point (per curve) after every five inte- 
gration steps. These values can be modified as indicated in the 
previous section. A run is terminated if the run exceeds 450 
lines or if more than 900 graph points are generated per curve. 

The accuracy of the numerical integration is controlled by 
the choice of step size. 'Unfortunately there is no simple 
method for making this choice other than by trial and error. 

As a first guess, a step size of 1/1000 of the total range of 
interest of the independent variable 'T' is usually reasonable. 
If one or more of the variables change rapidly in part of the 
range of 'T' , then a finer step size may be necessary for all 
or part of the run. The only way to be reasonably sure of the 
solution accuracy is to re-run the problem with a smaller step 
size and, then, to compare the results. A maximum of 4500 
integration steps is permitted in any one run. 

Data cards are assembled as follows: 

First: The user's job identification label is 



Second: 



punched in columns 1-42. 

The number of runs to be processed, which 
must be less than ten, is punched in column 



one. 



The run nimber together with the user's 



job identification label is placed on all 



output. 



Third: 



The order of the differential equations. 



which must be less than 31, is punched 



with two digits in columns one and two - 
for example, '03'. 
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Fourth: 



I 

I 



Fifth: 



Sixth: 



Seventh: 



Initial and final values of the independent 
variable (TI and TF) and the step size (DT) 
appear in the order TI - DT - TF and are 
punched, with decimal points, in columns 
1-10, 11-20 and 21-30. Also, it is possible 
to process the integration in either two or 
three segments of differing step size. The 
data, in columns of ten as above, then takes 
the form TI - DTI - TFl - DT2 - TF or TI - 
DTI - TFl - DT2 - TF2 -DT3 - TF . 

The values of the constants C(l), C(2), 

C(8) are punched, with decimal points, in 
columns 1-10, 11-20 ,..., 71-80 . Ten blank 
columns may be used for a zero or unused 
constant. If no constants are used, this 

entire card will be blank but it must 

not be omitted. 

The initial values of X(l) , X ( 2) , . . . , X (N) 
are punched, with decimal points, in 
columns 1-10,11-20,... . Additional cards 
are required if 'N' is greater than eight. 

Ten blank columns may be used for a zero 
initial value. 

This card controls the choice of variables 
for printout. Each group of ten columns, 

1-10 , 11-20 ,.. .up to 71-80, may be used to 
specify a column heading (eight characters) 
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Eighth: 



and a two-digit, right-justif ied subscript 
identifying the corresponding variable - 
for example: TIME 00 and SPEED 03 in 

columns 1-10 and 11-20 would cause the 
independent variable 'T', represented on 
this card by the subscript 00, to be 
printed out under the column heading 'TIME' 
and the variable X(3) to be printed out 
under the heading 'SPEED'. No more than 
eight variables can be output during one 
run. If no printout is desired, this card 
must be blank. 

This card controls the choice of variables 
for graph output. Up to four curves can be 
plotted, either all on separate graphs or 
all on a single graph. Each group of 20 
columns 1-20, 21-40 ,... specifies a curve - 
for example: SPEED VS. TIME 0300 (in columns 

1-20) would cause X(3) to be plotted 
vertically against the independent variable 
'T' horizontally. (Again, the subscript 00 
represents the independent variable 'T'.) 

The graph output would be labeled ' SPEED VS . 
TIME'. Note that the first 16 columns of each 
group are used for the label, the seventeenth 
and eighteenth for the 'Y' ordinates, and the 
last two for the ' X ' ordinates . 



66 



If the curves are to be on separate graphs , each graph 
must have a label. If the curves are to be all on one graph, 
only the first label, columns 1-16, must be provided. The 
other labels must be blank. The entire card must be blank if 
no graphical output is required. 

If several solutions are required of the same equations 
but with different constants or initial conditions, the number 
of runs, specified on the second data card, can be increased. 

For all runs after the first, only the fourth through the last 
data cards are to be supplied. Except for the data on the 
first three cards , and the equations themselves , no information 
is retained between runs. 

Caution should be exercised by the user because too large 
an integration step size can result in unacceptable errors in 
the solution, and even in instability in the integration process. 

This subroutine is maintained for use with the Naval 
Postgraduate School, Monterey, California, IBM 360 computer. 

It was programmed by J. R. Ward in October, 1963; was revised 
in November, 1971; and was redocumented in October, 1973. 
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